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ABSTRACT 

Globular star clusters are among the first stellar populations to have formed in the Milky Way, 
and thus only a small sliver of their initial spectrum of stellar types are still burning hydrogen on 
the main-sequence today. Almost all of the stars born with more mass than 0.8 Mq have evolved 
to form the white dwarf cooling sequence of these systems, and the distribution and properties of 
these remnants uniquely holds clues related to the nature of the now evolved progenitor stars. With 
ultra-deep HST imaging observations, rich white dwarf populations of four nearby Milky Way globular 
clusters have recently been uncovered, and are found to extend an impressive 5-8 magnitudes in the 
faint-blue region of the H-R diagram. In this paper, we characterize the properties of these population 
II remnants by presenting the first direct mass measurements of individual white dwarfs near the tip 
of the cooling sequence in the nearest of the Milky Way globulars, M4. Based on Gemini/GMOS and 
Keck/LRIS multiobject spectroscopic observations, our results indicate that 0.8 Mq population II 
main-sequence stars evolving today form 0.53 ± 0.01 Mq white dwarfs. We discuss the implications 
of this result as it relates to our understanding of stellar structure and evolution of population II stars 
and for the age of the Galactic halo, as measured with white dwarf cooling theory. 
Subject headings: globular clusters: individual (M4) - stars: evolution - stars: mass loss - techniques: 
spectroscopic - white dwarfs 



1. INTRODUCTION 

The end products of stellar evolution for 98% of all 
stars will be white dwarfs. This final state results be- 
cause most stars are not massive enough to ignite C and 
O, so that nuclear reactions cease with the formation of 
a degenerate core of He or a combination of He, C, and 
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O. With the cessation of nuclear reactions, white dwarfs 
emit light only via their cooling and become dimmer as 
time passes. These remnants have very thin surface lay- 
ers of H and/or He, the bulk of the envelope having been 
blown away during relatively quiescent stages of mass loss 
in the progenitor star. For the oldest stellar populations, 
>15% of the entire present day mass of the system is 
tied up in white dwarfs, and therefore their distribution 
and detailed properties represent a unique link to explore 
the nature of the now-evolved population II progenitors. 
Such studies are enhanced in star clusters where all of 
the remnants are co-spatial, have the same total age, and 
evolved from stars with the same metallicity, yet over a 
range in initial mass. 

Milky Way globular star clusters have been dated to 
a formation time roughly 12 Gyr ago (e.g., Krauss & 
Chaboyer 2003), and therefore represent the first struc- 
tures to form in the Galaxy. These systems are one of 
the most important tools that we have to probe the for- 
mation and evolution processes of galaxies (Brodie & 
Strader 2006). Given their ages, the present day main- 
sequence turnoff mass in globulars is about 0.8 Mq, 
and therefore the bulk of their population has evolved 
through post main-sequence evolution, the majority hav- 
ing formed faint white dwarf stars. Through almost 400 
orbits of Hubble Space Telescope imaging observations 
this decade, our team has uncovered the complete white 
dwarf cooling sequences of the two nearest globular clus- 
ters, M4 and NGC 6397 (Richer et al. 2004, 2008), and 
will target the third nearest cluster, 47 Tuc, in Cycle 17 
(GO-11677). A shallower study of the massive globular 
Omega Cen by Monelli et al. (2005) and Calamida et al. 
(2008) has successfully uncovered over 6500 remnants in 
this single cluster^" , roughly a quarter of the total known 

Given the larger distance, the brightest faint-blue objects on 



2 



Kalirai et al. 




12 



16 



20 



24 



28 



' I " " I " " I " " I ' 
NGC 6397 



lIi 



1 2 3 
V-I 




14 



16 



18 - 



20 - 



22 - 



24 



26 



28 



Omega Cen 




-10 12 3 4 
B-R 



Fig. 1. — HST CMDs of four nearby globular star clusters, extending from bright giant stars to the faintest white dwarfs (M4 - Richer 
et al. 2004, NGC 6397 - Richer et al. 2008, 47 Tuc - J. Anderson, 2009, private communication. Omega Cen - Calamida et al. 2008). 
The faint-blue remnants in each cluster represent evolution of stars that were more massive than the present day turnoff, M ~ 0.8 Mq. 
As discussed in § 1, the distribution and properties of these white dwarfs hold important clues for our understanding of stellar evolution 
processes in population II systems. 

the mass of the brightest white dwarfs that are forming 
today in these old systems. Yet, at present, no direct 
measurements of the masses of such population II white 
dwarfs exists, although several indirect arguments sug- 
gest that the value should be Mfinai ~ 0.53 Mq. As 
discussed in Renzini & Fusi Pecci (1988) and Renzini 
et al. (1996); stellar evolution theory and our under- 
standing of the maximum luminosity of red giant (RG) 
and asymptotic giant branch (AGB) stars, as well as the 
luminosity of the horizontal branch (HB), predict that 
white dwarfs forming in globular clusters today should 
have 0.51 < Affinal < 0.55 Mq. The best observational 
verification of this theoretical prediction in the case of 
cluster white dwarfs was provided recently by Moehler 
et al. (2004), who reported the masses of white dwarfs in 
NGC 6752 to be 0.53 Mq. However, their spectroscopic 
measurements were only of sufficient quality to yield the 
temperatures of the white dwarfs, and therefore the pho- 
tometry from the CMD was leveraged to derive masses. 
Recently, a study of the masses of field white dwarfs was 
also presented from the SPY project (ESO SN la Progen- 
itor surveY). Pauli et al. (2006) find several hundred thin 
disk and thick disk white dwarfs, and seven candidate 



white dwarf population in the entire Milky Way from all 
other studies combined (e.g., Eisenstein et al. 2006, and 
subsequent Sloan Digital Sky Survey work). The color- 
magnitude diagrams (CMDs) of these clusters, extending 
to impressively faint magnitudes, are presented in Fig- 
ure 1 (see also Bedin et al. 2009 for a more recent ACS 
study of M4). 

The white dwarf cooling sequences of nearby globu- 
lar clusters have been analyzed to independently (e.g., 
of main-sequence stars and turnoff morphology) derive 
the distances and ages of the clusters (Hansen et al. 
2004, 2007), explore age variations between metal- rich 
and metal-poor systems using stars whose evolution is 
almost completely independent of metallicity (Hansen 
et al. 2007), possibly suggest large fractions of helium- 
core white dwarfs produced through binary interactions 
(MoncUi et al. 2005; Calamida et al. 2008), and hint that 
mass loss stages leading to the formation of white dwarfs 
may involve small asymmetric kicks (Davis et al. 2008). 
For each of these studies, the interpretation of the pho- 
tometry to yield the results depends on knowledge of 

the apparent Omega Cen white dwarf cooling sequence may be 
subdwarf stars. 
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halo objects based on kinematic analysis. The masses 
of six of these white dwarfs range from 0.44 - 0.51 Mq, 
with an additional lower mass star at 0.35 Mq. 

In this paper we present the first results from a new 
study aimed at directly (i.e., purely from spectroscopic 
absorption line fitting) measuring the masses of known 
population II white dwarfs in different environments. 
Our goal is to target the brightest white dwarfs in the 
nearest globular star clusters (such as M4, NGC 6397, 
and 47 Tuc) and in the halo field, with ground-based 
spectroscopic instruments on 8-10 m telescopes. The 
clusters listed above span an appreciable range in metal- 
licity (factor of 20) and are located at distances such 
that the brightest white dwarfs with My = 10.5 have 
apparent magnitudes of y = 22 - 24. The faintness of 
these stars, crowding issues in the cluster, and the need 
to obtain high signal-to-noise (S/N) spectra of the higher 
order Balmer lines at A < 4000 A make this a difficult, 
but possible, project with current telescopes. Below, we 
discuss the results from our first study, in which we target 
two dozen white dwarfs in M4 with the Gemini/GMOS 
and Keck I/LRIS multiobject spectrographs over several 
observing runs. The spectra for six of these stars, that 
are confirmed as cluster members and that have suffi- 
cient S/N to yield accurate mass measurements, are con- 
sistent with Affinal = 0.53 ± 0.01 M©. The impact of this 
measurement for bettering our knowledge of stellar struc- 
ture and evolution, and for dating halo globular clusters 
through white dwarf cooling theory, are discussed. An 
analysis of the remaining M4 spectra that did not yield 
masses, e.g., to define the spectral types of population 
II white dwarfs, is presented in the companion paper by 
Davis et al. (2009). 

2. IMAGING AND SPECTROSCOPIC 
OBSERVATIONS OF M4 

The deep HST imaging data of M4 illustrated in Fig- 
ure 1 (left) were collected with the Wide Field Planetary 
Camera 2 (WFPC2) camera in 2001 as a part of GO-8679 
(Richer et al. 2004). The observations spanned a total 
of 123 orbits, in the broadband F606W^ and F^lAW fil- 
ters, and were centered on a field located at (q!,i5) = 
(16'^23'"54.6% -26°32'24.3"), 2.4 pc (4.8') S-E of the 
center of the cluster. Although the HST photometry is 
very accurate, and there is no ambiguity on the mem- 
bership of the white dwarfs, the bulk of the stars are too 
faint to be useful for spectroscopic follow up with mod- 
ern day facilities. Additionally, the WFPC2 field of view 
is 5.3 square arcminutes, about six times smaller than 
our ground based instruments, and therefore the spectro- 
scopic program would be highly inefficient with just this 
input list of targets. Below we first discuss supplemental 
imaging of these data with the Gemini South telescope 
(§2.1), followed by spectroscopic follow up (§2.2). In 
§§ 2.3 and 2.4, we discuss additional spectroscopic obser- 
vations of M4's white dwarf population with the Keck I 
telescope. 

2.1. Gemini/GMOS Photometry and Astrometry 

We reobserved M4 with the Gemini South telescope 
and Gemini MultiObject Spectrograph (GMOS), anchor- 
ing the deep HST field near the south-eastern edge of the 
ground based image and placing shallower HST point- 
ings near the western side. The GMOS photometry, over 
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Fig. 2.— The footprint of the HST/WFPC2 (dotted line), Gem- 
ini/GMOS (solid line), and Keck/LRIS (dashed) instruments. The 
Cycle 9 deep HST observations shown in Figure 1 comprise the 
bottom-right pointing, with the two other HST pointings being 
archive observations that we retrieved. The white dwarf candi- 
dates selected from the GMOS photometry (i.e., those stars with 
magnitudes between V = 20 and 24, and colors less than V — I ^ 
0.5) are shown as dots. The objects targeted for GMOS spec- 
troscopy are shown as diamonds and those for LRIS spectroscopy 
are represented with pluses. The center of M4 is located at (a, 5) 
= (254.898° ,—26.526) and is indicated with an asterisk. 

a 5.5' X 5.5' area, has the added advantage that it ex- 
tends continuous imaging into the denser region towards 
the center of M4 and therefore provides a higher density 
of targets down to a given magnitude limit. Of course, 
crowding in the ground-based data prevents the target- 
ing of white dwarfs near bright stars, and this increases 
in the direction of the cluster center. A footprint of the 
HST and GMOS images is provided in Figure 2. 

As this map illustrates, less than half of the Gem- 
ini/GMOS field of view overlaps the HST/WFPC2 ob- 
servations. In order to effectively select white dwarf 
candidates over the entire GMOS field, we obtained 
imaging in both the g' and r' filters, down to a depth 
of ^25th mag. The images were corrected for bias 
and flat fielding by the Gemini iraf pipeline, version 
1.4. The processed images were then reduced with the 
standard DAOPHOT/ALLSTAR photometry package to 
yield PSF photometry (Stetson 1994). Further details 
on the steps involved in our application are provided in 
Kalirai et al. (2008) for a similar study. 

The matched catalog of WFPC2 and GMOS sources 
contains 10'^ stars. We calibrate the g' and r' magni- 
tudes to HST F555W and F8UW (hereafter called V 
and /) using derived transformations for the common 
stars in the space- and ground-based systems (see Davis 
et al. 2009 for details on the transformations). These off- 
sets are then applied to the entire GMOS catalog. The 
positions of all stars on the GMOS image were carefully 
measured with respect to both HST positions and stan- 
dard USNO guide stars to yield precise astrometry. 

The Gemini/GMOS CMD of M4 is presented in Fig- 
ure 3. Both the main-sequence and white dwarf cooling 
sequence are well measured, the latter containing about 



4 



Kalirai et al. 



16 



18 - 



20 



22 



24 




LRIS 9 (CMOS 14) 



CMOS photomptry 



GMOS WD candidates 



0.5 1 



1.5 
Y-I 



2.5 



Fig. 3.— The CMD constructed from the CMOS photometry 
(small dots). There is a clearly identified white dwarf cooling se- 
quence extending from 22 < V < 24.5 with an approximate color 
of y — 7 = 0.3 (one brighter object is also present). This CMD is 
contaminated by a more distant field population, located blueward 
of the main sequence. There is also contamination in the white 
dwarf region of the CMD from field white dwarfs and possibly blue 
compact galaxies, which we address in § 4.3. The open circles en- 
circling the dots represent those objects which were targeted by 
Keck/LRIS for spectroscopy (see §2.4). The red open circles are 
additional objects that were targeted by Gemini/GMOS, and that 
are not in the Keck sample due to slit conflicts (see §§ 2.2 and 2.4). 
The open circles lacking small dots are HST detections for which 
we do not have ground based photometry (two of these stars are 
off the plot). These latter objects were selected for the sole pur- 
pose of flUing the mask in relatively empty regions, and therefore 
no pre-sclection was used. 

two dozen candidates brighter than V — 24. The loca- 
tion of these white dwarf candidates are illustrated on 
Figure 2 as smaller points. The selection here includes 
stars with 20 < V < 24.5 and colors V - I < 0.5. Note 
the paucity of ground based candidates in the two inner 
WFPC2 fields (top-left). There are in fact many white 
dwarfs at these positions, however, the crowding limits 
their detection from the ground. 

We further subdivide the sample of white dwarf can- 
didates into two bins, with priority 1 stars being our 
strongest white dwarf candidates in the dominant part of 
the cooling sequence, and priority 2 objects being some- 
what more uncertain objects that are located off the cool- 
ing sequence (e.g., redder or bluer) and/or near a bright 
star. These latter objects are only included in the spec- 
troscopic mask if they do not overlap a priority 1 star 
(hence, we get a spectrum of the potential white dwarf 
for free). Any of the confirmed members of M4, based 
on HST proper motions in the region of overlap, were 
also targeted as priority 1 objects. Finally, we include 
several redder stars that are candidate CVs in the pri- 
ority 2 sample; most of these are HST proper motion 
members as well. In total, 21 white dwarf candidates 
were targeted with Gemini/GMOS for follow up spec- 
troscopy. Note, a few of the very red or faint objects 
shown in the CMD were included without pre-selection 
for the purpose of filling the spectroscopic mask (these 




Fig. 4.— The CMOS photometry (upper left), WFPC2 pho- 
tometry (lower left), astrometry showing position relative to the 
instrumental footprints (upper right), and postage stamp section 
of the CMOS image centered on the target (lower right) for object 
WD 09. 

objects were not detected from the ground). 

An example of the GMOS and WFPC2 photometry, as- 
trometry, and a postage stamp image of the white dwarf 
candidate WD 09 is shown in Figure 4. Similar figures 
for all other white dwarf candidates are provided in the 
Appendix. 

2.2. Gemini/GMOS Spectroscopy 

In total, we obtained approximately 14.5 and 9 
hours of queue scheduled spectroscopic science expo- 
sures with Gemini/GMOS in June/ August 2005 and Au- 
gust/September 2006, respectively. The observations 
were taken with the B600 grating, covering a broad wave- 
length region of ~2700 A centered at 4680 A. The data 
were binned by a factor of two in the spectral direction to 
improve the S/N. All science exposures were obtained in 
sub-arcsecond seeing conditions and photometric skies. 

The raw data frames were downloaded from the Cana- 
dian Astronomy Data Center (CADC) in multi-extension 
FITS (MEF) format and reduced with the Gemini iraf 
Package, version 1.4. The reductions are described in 
detail in Davis et al. (2009). Of the 21 targeted white 
dwarfs, a spectrum was recovered for all but three. Re- 
markably, all 18 of the recovered spectra show pressure- 
broadened absorption lines indicative of white dwarfs. 
The spectra for these 18 stars, as well as an analysis 
of their spectral types is also presented in Davis et al. 
(2009). Unfortunately, in order to determine reliable 
spectroscopic masses, characterization of He is crucial, 
and even H8 (3889 A) is very useful. The lack of flux at 
these wavelengths in our Gemini data limited the utility 
of the GMOS spectra for the central goals of this paper. 
Although we could derive mass measurements from the 
spectra (see §4.1), the uncertainties are several tenths 
of a solar mass and therefore we decided to supplement 
these data with higher S/N observations, especially at 
bluer wavelengths. 
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2.3. Keck/LRIS Target Selection 

With a sample of confirmed white dwarfs established 

from the GMOS observations, we apphed for Keck I tele- 
scope time to rcobscrvc the targets using the LRIS multi- 
object spectrograph. Not only docs Keck provide a larger 
aperture than Gemini, but the spectrograph is very effi- 
cient in the blue and has an atmospheric dispersion cor- 
rector to minimize light loss at large airmasses (this is 
important for this study since Keck is in the north and 
M4 is in the south). Given the similarities in the field 
of view (LRIS has a 5'x 7' footprint), we simply used 
as many of the confirmed sample of white dwarfs from 
Gemini as our input, as well as several additional redder 
stars in the same brightness range. The location of the 
31 LRIS targets in the ground-based CMD are shown in 
Figure 3 as larger open circles. The objects marked in 
red circles arc those that were targeted with Gemini, but 
not Keck. The LRIS footprint is also shown in Figure 2. 

2.4. Keck/LRIS Spectroscopy 

Keck/LRIS spectroscopic observations (R. M. Rich, 
PI) of the single M4 mask proved challenging. Half 
nights, required due to the low declination of M4, were 
scheduled in 2005 June, 2007 April/July, and 2008 April, 
with the 2007 and 2008 observations benefiting from the 
LRIS atmospheric dispersion correction system. The 
multiple allocations were required due to poor weather 
that resulted in the loss of a significant number of sched- 
uled hours. Altogether, we obtained 10.6 hours of useful 
exposure time. Even these data are of highly variable 
quality, with seeing measurements ranging from 0"8 to 
over 2'.'0. In these crowded fields, our final spectra are 
dominated by a few of the contributions from the best 
nights. 

LRIS (Oke et al. 1995) has a dichroic that splits 
the spectrum into two channels at -^5500 A. On the 
blue side, we use the 400/3400 grisni (dispersion = 1.09 
A/pixel) which covers a wavelength region from the at- 
mospheric cutoff to the dichroic. For the red side, we use 
the 600/7500 grating (dispersion = 1.28 A/pixel), cen- 
tered at 6600 A, which covers a wavelength baseline of 
2620 A. At our resolution and central wavelength, only 
Ha, which is a rather poor mass indicator, landed on the 
red side. Hence, we chose not to reduce the red-side data 
for this paper. 

The blue-side LRIS spectra were reduced using stan- 
dard IRAF tasks. The trace, sky subtraction, and ex- 
traction were all performed with the APALL task. As 
expected given the better blue sensitivity of LRIS com- 
pared to GMOS, the trace at blue wavelengths was far 
more certain. Given the high levels of crowding and back- 
ground light from bright M4 stars in the field, our abil- 
ity to obtain a reliable trace, and perform accurate sky 
subtraction is dependent on the particulars of each indi- 
vidual slitlet (see §3). 

The wavelength calibration was calculated in iRAF 
from spectra of three lamps containing Hg, Zn, and 
Cd. The typical final dispersion of the wavelength solu- 
tion was well constrained, with uncertainties of ~0.1 A. 
The reduced, wavelength calibrated LRIS spectra were 
flux calibrated using the spectrophotometric standard 
HZ44. Finally, the individual spectra were combined us- 
ing SCOMBINE, with weights assigned according to the 



TABLE 1 

Properties of the M4 White Dwarf Targets 



Identification # 
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These stars were not detected from the ground; magnitudes are 
from the HST photometry. 



individual S/N ratios. 

We summarize the identification, positions, and bright- 
nesses of all white dwarf candidates in Table 1. The four 
stars that have slit conflicts, and therefore were not ob- 
served with Keck, are listed at the bottom of the table. 
Altogether, 31 objects were targeted with Keck/LRIS. 
Note, the identiflcation of these stars as "WD" at this 
stage does not confirm their nature, although we will 
find later that most are in fact white dwarfs. We use 
this naming scheme for convenience as we will refer back 
to this Table in future sections. 

3. ELIMINATING NON- WHITE DWARFS AND STARS 
WITH POOR SPECTRA 

The ground-based CMD of our targets in Figure 3 

shows that the dominant, blue white dwarf cooling se- 
quence of M4 extends from V = 22.5 to 23.5, and con- 
tains eleven stars that were spectroscopically observed 
with LRIS (note, the clump at F ~ 22.7, V-I ~ 0.3 con- 
tains five LRIS targets). An additional target, WD 02, is 
almost one magnitude brighter than the family of fainter 
points and has V = 21.3, another four stars are found 
fainter and bluer than the cooling sequence at V ^ 23.5, 
y — J ~ 0.25 (including the one object without ground 
based photometry), and two stars are at the faint end 
and redder than the cooling sequence at V = 23.2 and 
V — I = 0.4. Most of the stars that are redder than this 
sample of 18 objects (i.e., V — I > 0.5) yielded spectra 
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Fig. 5. — Keck/LRIS spectra for ten white dwarfs along the sightline to M4 (left). The S/N per resolution element at a wavelength of 
~4030 A in these data is 30 - 50; the higher order hydrogen Balmer lines are well characterized. The spectra for each of these stars is fit to 
atmosphere models to yield T^g, log g, and Affinal for the white dwarfs as described in § 4.1 (top-right). In these smaller panels, the Balmer 
lines for a given white dwarf are displayed with H/3 at the bottom, and successively higher order transitions on top. The simultaneous fit 
to all Balmer lines in the best fit model is shown as a smooth red curve. A closer view of the white dwarf cooling sequence of M4 is shown 
in the CMD in the bottom-right panel, with a 0.5 Mq theoretical cooling model (Wood 1995). Note, all spectra in this figure have been 
smoothed by a three pixel boxcar function for clarity. 



that are not interesting for the goals of this paper and 
will be commented upon further in Davis et al. (2009). 
In Table 1, these objects are WD 01, 03, 07, 08, 10, 11, 
13, 14, 21, 23, 26, 27, and 30. The one exception to this 
is possibly the relatively bright star just beyond the V — I 
= 0.5 color cut atV = 22.54, V - I ^ 0.51 (WD 23). 
This star is a proper motion member of M4 but, as shown 
in the Appendix in Figure A3, is located very close to a 
bright neighbor. Likely, the color (and spectrum) of this 
star is contaminated from the neighboring star and we 
therefore eliminate it from our analysis. 

The family of 18 white dwarf candidates in M4's 
'pseudo' cooling sequence are all confirmed to be white 
dwarfs (see Davis et al. 2009). For our primary goal 
of measuring accurate masses, we want to first restrict 
this sample to only include well measured, isolated white 
dwarfs with robust S/N (see §4.1). The stars along 
the cooling sequence in our sample vary in brightness 
by about a factor of three, and therefore the fainter 
stars will intrinsically have much less signal than their 
brighter counterparts. Added to this, we note that the 



final quality of the spectrum for any star will depend on 
several factors and therefore not always correlate with 
the brightness. For example, the spatial length of the 
given slit effects the quality of the sky subtraction, the 
position of the star relative to M4's center effects the sky 
background level, and any astrometric mis-centering of 
the star within the slit causes light loss. We first elim- 
inate three (fainter) stars that are located within <30 
pixels (~2") of a bright neighbor. These stars, WD 16, 
19, and 22, all exhibit a noisy spectrum that is likely con- 
taminated with light from the nearby, brighter star (see 
Figures A2 & A3 in the Appendix). We also eliminate 
five white dwarfs, four of which are isolated, for which 
our extractions resulted in S/N < 10 spectra, WD 12, 
17, 18, 25, and 28. The extracted spectra for these white 
dwarfs suffer from unusable line profiles blueward of H^, 
but with improved S/N, these particular sources could 
also be analyzed as described below. 

Our final sample has ten white dwarfs, all of which 
are close to the M4 cooling sequence and have V < 23.4 
(seven of which have V < 23). As we demonstrate in 
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ID 


V 




log .9" 




T/ ^ 

•'theory 


^cool 




Member'' 


WD 00 


23.32 


20,600 ± 600 


7.75 ± 0.09 


0.52 ± 0.04 


22.82 ± 0.26 


40 ± 6 


39 


YES (PM) 


WD 04 


22.69 


24,300 ± 500 


7.68 ± 0.08 


0.50 ± 0.03 


22.40 ± 0.23 


19 ± 1 


47 


YES (Lum) 


WD 06 


22.65 


25,600 ± 500 


7.87 ± 0.08 


0.59 ± 0.04 


22.59 ± 0.23 


18 ± 2 


51 


YES (Lum) 


WD 15 


22.73 


24,300 ± 600 


7.79 ± 0.09 


0.55 ± 0.04 


22.57 ± 0.27 


21 ± 2 


29 


YES (Lum) 


WD 20 


23.01 


19,700 ± 600 


7.73 ± 0.10 


0.51 ± 0.05 


22.89 ± 0.29 


49 ± 8 


32 


YES (PM) 


WD 24 


22.72 


25,700 ± 500 


7.70 ± 0.07 


0.51 ± 0.03 


22.33 ± 0.21 


16 ± 1 


54 


YES (PM) 


WD 02 


21.32 


18,800 ± 200 


8.16 ± 0.04 


0.75 ± 0.03 


23.56 ± 0.14 


128 ± 11 


61 


NO 


WD 05 


22.71 


28,200 ± 400 


7.59 ± 0.08 


0.48 ± 0.03 


21.93 ± 0.22 


12 ± 1 


49 


NO 


WD 09 


22.50 


25,100 ± 600 


8.19 ± 0.09 


0.79 ± 0.05 


23.10 ± 0.25 


44 ± 12 


32 


NO or DD 


WD 29 


22.94 


20,900 ± 700 


7.44 ± 0.11 


0.41 ± 0.04 


22.32 ± 0.31 


33 ± 4 


27 


NO 



"The new Stark broadening calculations in Tremblay & Bergeron 
(2009) have not been used to recalculate the Teff and log g values 
(which would require approximate shifts of +1,500 K and 0.1 dex), 
although the masses have been corrected as discussed in §4.1. 

Theoretical luminosity from spectral fits (sec §4.3). The error 
bar includes the la error in the distance modulus (0.09 mags). 

The spectral S/N per resolution element, calculated between H5 
and He. 

Membership evaluation (see §4.3); PM = Proper Motion Mem- 
ber, Lum = Luminosity Member, DD = Double Degenerate Can- 
didate. 

§ 4.3, four of these ten white dwarfs, WD 00, 09, 20, 

and 24, are HST proper motion members of M4. We 
present the spectra for these ten white dwarfs in Figure 5 
(left) . All of the spectra clearly show pressure-broadened 
Balmer lines extending from H/3 to H8. These are there- 
fore all DA (hydrogen) spectral type white dwarfs. We 
report the spectral types of all of our white dwarfs, in- 
cluding stars that are too faint for mass constraints, in 
Davis ct al. (2009). The S/N per resolution clement in 
our ten stars is generally 40 - 50, and at best 60 for one 
star (WD 02). 

We also present a detailed view of the white dwarf 
cooling sequence in the CMD, and indic;atc the positions 
of these ten white dwarfs in the bottom-right panel of 
Figure 5 (dotted circles). The remaining open circles 
(without dots) in this diagram are those stars that were 
spectroscopically targeted in this narrow region of CMD 
space. We note that the magnitudes and colors of WD 15 
and 24 are almost identical and so the two dots can not 
be discerned. We also illustrate a theoretical white dwarf 
cooling model for M = 0.5 Mq (Wood 1995), which cuts 
through our sample nicely. From this plot, it is clear 
that WD 02 is well removed from the dominant sequence 
of points and the cooling model. Although it is possible 
that the photometry for this star is in gross error, it is 
more likely that the star is a field white dwarf along M4's 
sightline (see §4.3 for more information). 

4. MASS MEASUREMENTS OF WHITE DWARFS IN M4 

The technique of measuring masses of white dwarfs 
from Balmer line fitting has been widely employed for 
field stars, such as in the Palomar-Green (PG) Sample 
(Liebert, Bergeron, & Holberg 2005) and Sloan Digital 
Sky Survey (SDSS) (Eisenstein et al. 2006; Kepler et al. 
2007). The technique has also been cross-checked against 
independent measures of the mass of a white dwarf such 

The spectral S/N is calculated at ~4030 A as discussed in 

§4.2. 



as astrometry for binary systems (e.g., Sirius B, Barstow 

ct al. 2005); fits to the mass-radius relation for stars with 
trigonometric parallaxes (Holberg, Bergeron, & Gianni- 
nas 2008), the gravitational redshift method (Reid 1996), 
and pulsation mode analysis (Kawaler 1991). For tem- 
peratures hotter than ~ 12, 000 K, these methods are in 
excellent agreement with the spectroscopic line fitting 
technique (Bergeron, Liebert, & Fulbright 1995). 

In this section, we first briefiy describe our method 
to calculate masses for the ten white dwarfs discussed 
above. We then verify the uncertainties in the masses 
using synthetic spectra of the same qiiality as the obser- 
vations, and fit using the same methods described below. 
Finally, we remove field white dwarfs from the sample 
and analyze the mass of the confirmed M4 members. 

4.1. Balmer Line Fits 

Fitting the Balmer lines of a white dwarf to model 
atmospheres involves reproducing multiple line profiles 
from H/? to higher order transitions (Bergeron, Saffer, 
& Liebert 1992). Typically, four or five well measured 
lines are enough to constrain the temperature (Toff) and 
surface gravity (log g) of the star (Bergeron, Saffer, & 
Liebert 1992). The fifth Balmer line in this sequence, 
H8, is the weakest and bluest (at A = 3889 A), and there- 
fore a spectrograph with high throughput in the blue 
is needed to characterize the feature accurately in faint 
stars. General details of the fitting technique used to 
derive Toff and log g are provided in Bergeron, SaS^er, 
& Liebert (1992), and recent refinements to this method 
are given in Liebert, Bergeron, & Holberg (2005). Addi- 
tional information on the fitting of white dwarf spectra 
similar to these observations, and taken with our instru- 
mental setup, are presented in § 6 of Kalirai et al. (2008). 
Specifically, the improved method to normalize the line 
profiles described in Liebert, Bergeron, & Holberg (2005) 
is used in our work. 

Simultaneous fits of model atmospheres to all of the 
Balmer lines, in each of the ten white dwarfs discussed 
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above, arc presented in Figure 5 (top-right). Within each 
panel, the absorption hnes of a given white dwarf are 
arranged with H/3 at the bottom, and successive higher 
order transitions on top. The model fits are shown as 
smooth red curves and are excellent in all cases. We also 
confirmed that the two anomalous features, small spikes 
caused by bad sky subtraction in the H/3 line of WD 00 
and 29, do not affect the results of the fit. This line is 
least sensitive to changes in the gravity of the star. The 
output from the spectroscopic fits are the T^g and log g 
for each star, which we summarize in Table 2. 

The masses of these ten stars (Affinal) arc simply calcu- 
lated by combining the measured surface gravity of the 
stars with the well established mass-radius relation for 
white dwarfs. In general, most of the white dwarf masses 
are ~0.5 Mq with an uncertainty of <10%. We also es- 
timate both the cooling ages, tcooi, and the theoretical 
luminosity of each star by interpolating the Toff and log g 
within evolutionary models similar to Fontaine, Bras- 
sard, & Bergeron (2001), but assuming a 50/50 carbon- 
oxygen core mix.^^ Our default calculations assume thick 
hydrogen layers with q(R) = Mh/M ~ 10~^ and helium 
layers of (/(He) = 10^^, and are summarized in Table 2 
and discussed further below. We also tested the sensi- 
tivity of our results to models (e.g.. Wood 1995) with 
pure carbon cores and to models with very thin hydro- 
gen layers (e.g., (?(H) = 1G~^°). The results from these 
comparisons, for both the final remnant masses and the 
theoretical luminosities, are identical to our default val- 
ues within the measured uncertainties in these properties 
(which are small, see Table 2). We note, however, that 
the cooling ages of these bright, young white dwarfs are 
systematically larger in the pure C core models by 20 - 
30%. 

At the time of writing this paper, Tremblay & Bergeron 
(2009) have just presented the first results from an im- 
proved calculation of the Stark broadening of hydrogen 
lines in dense plasmas that are representative of white 
dwarf atmospheres. Their new results, including non- 
ideal effects, suggest that the masses measured through 
the spectroscopic technique require a slight upward revi- 
sion of 0.034 Mq. As the temperatures of the stars using 
the new calculations are also slightly hotter (e.g., 1,000 
- 2,000 K), the luminosities are unaffected. The full im- 
plications of this study have not yet been investigated 
and the dependency of this offset on mass is not well 
constrainc^d. In the results that follow, we have adjusted 
our default masses by +0.034 Mq to be consistent with 
these new calculations (including the reported masses in 
Table 2). We have not corrected the individual Tefi and 
log g values by a common offset (which would be approx- 
imately 1,500 K and 0.1 dex, respectively). 

4.2. Testing the Uncertainties in Mfinai' Synthetic 

Spectral Fits 

The S/N of the ten white dwarfs in our sample is listed 
in column 8 of Table 2. This is calculated from the 
pseudo-continuum between the B.6 and He Balmer lines, 
at ~4030 A. The pixel-to-pixel RMS scatter is converted 

to a S/N per resolution element by multiplying by the 
square root of the number of pixels per FWHM (six, in 

The models are available at 

http://www.astro.umontreal.ca/ bergeron/CoolingModels/. 



our case). The bulk of the stars have S/N = 40 - 50, 
although a few are as low as ^30 and one object is above 
S/N = 60. Note, these S/N values should only be used 
as a relative ranking of the quality of each spectra, and 
are not indicative of the global S/N. We can see that, 
due to the throughput, the S/N decreases blueward of 
the wavelength at which we made this measurement. 

In Figure 6 (top), we produce simulated white dwarf 
spectra over a range of S/N. These spectra have been 
convolved with a Gaussian profile with 6 A FWHM. Di- 
rectly comparing the quality of our data in Figure 5 to 
these simulated spectra, we can conclude that most of 
our objects are better than the S/N = 30 case, and simi- 
lar, if not better, than the S/N = 40 case. For the latter, 
we generate a Monte-Carlo with 500 trials and fit each of 
these spectra in our fitting code. The prescription used 
in the fitting is identical to that for our actual white 
dwarfs, namely, the same models are used in the fit with 
five Balmer lines. The results from this test are shown in 
the middle panel of Figure 6. The input stellar param- 
eters for the synthetic white dwarf are Toff = 25,000 K, 
log g = 7.80, and M ^ 0.536 Mq. From the distribu- 
tion of recovered masses, we find (M) = 0.54 Mq with a 
dispersion of 0.06 Mq. Therefore, we expect to measure 
the masses of our white dwarfs to roughly 10% precision 
based on this test, which agrees nicely with our reported 
uncertainties in mass (see Table 2, column 5). In the 
bottom panel of Figure 6, we extend the test above for 
other S/N values and plot the associated error from this 
fitting. 

4.3. Establishing Cluster Membership 

Thus far, we have established masses for ten white 
dwarfs in the M4 direction. As we have reported in previ- 
ous spectroscopic studies of white dwarfs in star clusters, 
some fraction of the stars along these lines of sight are ac- 
tually members of the field population and therefore need 
to be removed before discussing the mass distribution of 
bonafide members. The best method to establish such 
membership is through proper motion analysis, and we 
have such data from the HST overlap of a portion of our 
LRIS field (see Figure 2). Out of these ten white dwarfs, 
we can confirm that WD 00, 09, 20, and 24 are proper 
motion members (as are WD 16 and 23; these are re- 
jected from our sample based on the poor quality of the 
recovered spectra and/or contamination from a nearby 
neighbor - see §3). This is not a surprising result since 
the proper motions were folded into the selection process 
to target white dwarfs in the region of overlap with our 
HST field. 

In addition to these five confirmed M4 white dwarfs, 

we note that the remaining objects in our sample closely 
follow the 0.5 A'/q cooling sequence, and are consistent in 
position in the CMD with the confirmed members. Al- 
though we do not have proper motions for these stars, 
some of them are likely members of M4. To distinguish 
members from non-members, we can use the theoretical 
luminosity of each star to calculate its distance modu- 
lus. Any stars that have significantly discrepant values 
from the known distance of M4, (m — M)v = 12.51 ± 
0.09 (Richer et al. 2004), should be eliminated from our 
sample as they likely represent field white dwarfs along 
the line of sight, binary white dwarfs, or possibly objects 
for which the spectroscopic fit is incorrect due to some 
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Fig. 6. — Synthetic white dwarf spectra arc created over a range 
of S/N from 10 to 150, as measured at ~4030 A between the US 
and He Balmer hnes (top panel). The typical white dwarf in our 
observations shown in Figure 5 is similar in spectral quality to the 
S/N = 40 case in this plot. In the middle panel, wc illustrate the 
results from re-fitting 500 simulated spectra with this S/N. The 
input spectra in this Monte Carlo trial all had T^g = 25,000 K, 
log g = 7.80, and M = 0.536 Mq. Our results, using the same 
methods of fitting five Balmer lines as outlined in §4.1, indicates 
that the mean recovered mass is (M) = 0.54 Mq with a l-cr error in 
the mean of 0.06 Mq. In the bottom panel, we illustrate the error 
in the recovered mass from similar trials of 500 synthetic spectra 
for a series of S/N ratios. For values comparable to our data, we 
should be able to recover the masses of these white dwarfs to ~10%. 



unknown artifact in the spectrum. 

We compare the calculated distance modulus of each of 
our ten white dwarfs with the known M4 distance in Fig- 
ure 7. The 1(7 uncertainty in the theoretical magnitude 
(i.e., from spectral fitting) is plotted as a horizontal bar 
for each star, and the mean distance of M4 and ±2a un- 
certainties are indicated with the solid and dashed lines. 
We can use the proper motion members of the cluster as 
a rough guide to establish the expected scatter of known 
members in this diagnostic. The calculated distances of 
three of the four proper motion members of the cluster, 
WD 00, 20, and 24, are close to, but slightly larger than 
the distance of M4, whereas the fourth star (WD 09) ap- 
pears to be 0.7 mags closer than the cluster. This object 
could therefore be a non-member that shares the proper 
motion of the cluster. Alternatively, if the object is a 
member, its overluminous nature may suggest an unre- 
solved double degenerate in the cluster. We note that 
the mass derived from the spectrum of this star yields 
the unexpectedly high value of M = 0.79 Mq (see Ta- 
ble 2). Three additional white dwarfs are clearly mem- 
bers of the cluster, WD 04, 06, and 15, whereas three 
other stars WD 02, 05, and 29 are possibly field white 
dwarfs. 

To summarize, our final secure sample; of singly 
evolved, M4 member white dwarfs includes three objects 
that are proper motion members (WD 00, 20, and 24) 
and three objects which are on the cooling sequence and 
have the correct distance from spectral fitting (WD 04, 
06, and 15). 
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Fig. 7. — The combination of the theoretical magnitude of each 
white dwarf, from fitting the spectrum with models, and the ob- 
served magnitude is used to calculate individual distance moduli. 
The length of the bar for each white dwarf represents the la error 
in the theoretical magnitude. The solid line indicates the distance 
modulus of M4 (Richer et al. 2004) of (m - M)v = 12.51 ± 0.09, 
and the dashed lines represent 2cr bounds. The three objects at 
the bottom of the diagram, WD 00, 20, and 24, arc known proper 
motion members of M4. A further proper motion member, WD 09, 
is overluminous and appears to be 0.7 mags closer than the clus- 
ter. Likely, this star is either a field object sharing the proper 
motion of the cluster or an unresolved double degenerate. Another 
three objects in the middle of the diagram, WD 04, 06, and 15 are 
all consistent with membership whereas three objects at the top, 
WD 02, 05, and 29 are excluded from our member sample. 

4.4. Analysis o/Mfinai 

Theory predicts that the masses of white dwarfs form- 
ing in old, population II systems such as globular clusters 
are expected to be in the range 0.51 < M < 0.55 Mq 
(Renzini & Fusi Pecci 1988; Renzini et al. 1996). As Ta- 
ble 2 indicates, the six single white dwarfs that we have 
confirmed as members of M4 have masses that range from 
0.50 - 0.59 Mq. The uncertainties in our mass measure- 
ments are all roughly the same, under 10%, and therefore 
the un-weighted mean mass of the sample of six white 
dwarfs is measured to be Mfinai = 0.53 ± 0.01 Mq. 

The present day temperature of these cooling remnants 
is also well measured from our data. There is a group 
of four stars in our sample with V ^ 22.6 (My — 10.1) 
and all have the same temperature, T^e ~ = 25,000 K. 
The two fainter stars, WD 00 and 20, are cooler than 
this group by about 5,000 K, as expected. The surface 
gravities, log g, of all six stars are in the range log g = 
7.68 - 7.95. As mentioned earlier, the exact T^g and log g 
values have not been adjusted to reflect the small correc- 
tions required by the new Stark broadening calculations 
(Tremblay & Bergeron 2009). 

Our sampling of M4 cluster member white dwarfs in 
this study puts us in the fortunate situation of having 
characterized a group of four remnants at the same posi- 
tion on the CMD. Hansen et al. (2007) demonstrate that 
the mass difference between white dwarfs at the tip of 
the cooling sequence to those five magnitudes fainter at 
the hmit of the HST study is only ~0.1 Mq. These four 
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Fig. 8. — A S/N 100 composite spectrum of four individual 
M4 white dwarf members, all of which have the same luminosity 
to within 0.07 mags. The variation in mass between the individual 
stars is expected to be negligible. The spectroscopic fit to the 
Balmer lines of the composite spectrum yields a mass of M = 0.54 
± 0.02 Mq, in excellent agreement with our overall mean mass 
of the total confirmed sample of member white dwarfs. As for 
the spectra shown in Figure 5, we have smoothed the composite 
spectrum by three pixels in this Figure. 



white dwarfs, WD 04, 06, 15, and 24, differ in luminosity 
by just 0.07 mag, and therefore should have the same 
mass, and have evolved from stars with the same pro- 
genitor mass. The expected temperature range in this 
group of four stars should only be ^2%. As these dif- 
ferences are much smaller than our uncertainties, we can 
verify our mean mass measurement of M4 white dwarfs 
by creating a composite white dwarf spectrum, with very 
high S/N, from the coaddition of these individual spec- 
tra. Liebert et al. (1991) have shown that for unresolved 
double degenerate systems, the atmospheric parameters 
derived from assuming a single star are in fact an aver- 
age of the parameters of both components of the system. 
Therefore, spectra with identical gravities should yield 
the same average gravity in a coaddition as compared to 
the true individual values. 

Such a coaddition has the advantage that the resulting 
fit will be minimally influenced by any artifact affecting 
a given Balmer line in one spectrum (e.g., presence of 
a cosmic ray). This method has several disadvantages 
as well, for example, the spectral lines in the composite 
may be artificially blurred if the wavelength solution is 
not well determined in each of the individual spectra. We 
stress that this test is only intended as a verification of 
our results and should not be attempted in samples where 
the mass variation between individual white dwarfs is 
non- negligible (e.g., similar studies in younger clusters). 

The composite spectrum from this analysis is shown in 
Figure 8 (top), and exhibits a S/N - 100 at A = -4030 A. 
The Balmer lines in this spectrum are very accurately 
characterized; the resulting best fit model to the five lines 
is shown in the lower panel. The spectroscopic mass 
from the fit is measured to be M = 0.54 ± 0.02 Mq, in 
excellent agreement with our findings above. 



Overall, from our spectroscopic study of white dwarfs 
in M4, we can conclude that the mass of population II 
remnants forming today is Mgnai = 0.53 Mq and the 
temperature of stars presently at the tip of the white 
dwarf cooling sequence is Toff = 25,000 K. 

5. STELLAR EVOLUTION THEORY 

5.1. The Initial-Final Mass Relation 

The mass distribution of white dwarfs in nearby open 
star clusters has been mapped since 1977 when Volker 
Weidemann compared theoretical models of mass loss 
(e.g., Fusi-Pecci & Renzini 1976) to the observed masses 
of a few remnants in the Hyades and Pleiades clusters 
(Weidemann 1977). Since this pioneering work, similar 
mass measurements have been made for over 100 stars in 
a dozen open star clusters. A review of the earlier work 
is provided in Weidemann (2000) and a compilation of 
some more recent results is presented in Fcrrario et al. 

(2005) , excluding very recent studies by Dobbie et al. 

(2006) ; WiUiams & Bolte (2007); Kahrai et al. (2007, 
2008); WiUiams, Bolte, & Koester (2009); Dobbie et al. 
(2009). The mass distribution of the remnants varies 
as a function of the turnoff mass of their parent clus- 
ter, suggesting that more massive progenitors produce 
more massive white dwarfs (see Salaris et al. 2009 for a 
detailed analysis of the present relation). New measure- 
ments by our team in three older star clusters with ages 
of a few Gyrs have recently permitted, for the first time, 
a purely empirical fit to the data without the need for a 
theoretical anchor at low masses (see below). 

With the present study, we can once again extend the 
initial-final mass relation to new bounds. As the initial 
mass function in most (all?) stellar populations is bot- 
tom heavy, the characterization of mass loss from these 
low mass stars affects our understanding of stellar evolu- 
tion for the bulk of all stars today. In Figure 9, we present 
all of the constraints to date at the low mass end of the 
relation, corrected by -t-0.034 Mq from our previous re- 
porting to refiect the new Stark broadening calculations 
(see §4.1). The data includes the two confirmed clus- 
ter members in the 1.4 Gyr cluster NGC 7789 (Kalirai 
et al. 2008), the two white dwarfs in the 2.5 Gyr cluster 
NGC 6819 (Kalirai et al. 2008), the single carbon-oxygen 
core white dwarf in the 8.5 Gyr cluster NGC 6791 (Kali- 
rai et al. 2007), and the six white dwarfs in M4 from this 
study. For the latter, we have fixed the progenitor mass 
to the expected mass at the turnoff in an old, metal-poor 
population, Minitiai — 0.80 Mq, with an adopted error of 
0.05 Mq.^^ The data points for individual white dwarfs 
in these clusters are shown as open circles and binned 
averages of the population in each cluster are illustrated 
with larger, filled circles. For older clusters, all of the 
white dwarfs at the top of the cooling sequence evolved 
from progenitors just above the present day turnoff, and 
therefore we expect a single mass at the cooling sequence 
in a given cluster. Of course, in younger clusters, our 
spectroscopic measurements can include stars along the 
entire cooling sequence of the cluster and therefore there 

The exact turnoff mass depends on the metallicity and age 
of the cluster, and the choice of theoretical model assumed. For 
example, for alpha-enhanced models with [Fe/H] = —1.1 and t = 
12.0 Gyr, the VandenBerg, Bergbusch, & Dowler (2006) models 
yield a turnoff mass of 0.82 Mq. 
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Fig. 9. — The low mass end of the initial-final mass relation 
showing white dwarf mass measurements and calculated progenitor 
masses in NGC 7789 and NGC 6819 (Kalirai et al. 2008), NGC 6791 
(Kalirai et al. 2007), and the present study for M4. All white dwarf 
mass measurements have been corrected by +0.034 Mq to reflect 
the new Stark broadening calculations presented in Tremblay & 
Bergeron (2009). As discussed in §5.2 and Kalirai et al. (2007), 
the single white dwarf in NGC 6791 is much more massive thaji the 
mean mass of remnants in the cluster, most of which are low mass 
helium-core white dwarfs (and not shown here). The individual 
white dwarf - progenitor pairs are shown as open circles and the 
binned averages in each cluster are shown as larger filled circles. 
The relation shows a linear decline of final remnant mass with 
decreasing initial progenitor mass, extending down to the extreme 
low mass end, anchored by our M4 measurement at Mjnitiai = 0.80 
± 0.05 Mq and Affinal = 0.5.3 ± 0.01 Mq. The solid line shows the 
best fit empirical relation to all individual data points excluding 
the new M4 white dwarfs, and the dashed curve shows the mass of 
the core at the first thermal pulse (sec §5.1). 

can be a largo dispersion in white dwarf masses in a 
single cluster (e.g., NGC 2099, Kalirai et al. 2005 and 
NGC 2168, Williams, Bolte, & Koester 2009). 

The low mass end of the initial-final mass relation in 
Figure 9 now contains a total of eleven data points, and 
shows a roughly linear decrease in white dwarf mass with 
initial mass, anchored at the extreme left end by the new 
M4 measurement, Minitiai = 0.80 ± 0.05 Mq and Mgnai = 
0.53 ± 0.01 Mq. Given that the oldest observable stellar 
populations in the Universe are the same age as M4, this 
relation is now complete at the lower limit and can not 
be observationally extended further to lower masses. As 
discussed earlier, several predictions of stellar evolution 
theory have long suggested that the masses of population 
II white dwarfs forming today should be 0.51 < Mfinai < 
0.55 Mq (Renzini & Fusi Pecci 1988; Renzini et al. 1996), 
and therefore our first direct measurement is in perfect 
agreement with these predictions. Based on our mea- 
surements, the total integrated post main-sequence mass 
loss of low mass stars in old population II systems will 
be -34%. 

Theoretically, there is an expected lower limit to the 
masses of carbon-oxygen core white dwarfs. This is given 
by the mass of the stellar core at the first thermal pulse 
on the asymptotic giant branch, an evolutionary point 
before the core has had a chance to grow on the branch. 
This core mass has a very small dependency on the stellar 



metallicity. In Figure 9, we illustrate this expected core 
mass as a dashed curve (Girardi et al. 2000). Over the 
range of initial masses on this plot, the core mass is ex- 
pected to remain roughly the same at Mgnai = 0.52 Mq, 
showing a flat relation for Minitiai ^ 1 - 2 Mq. At 
0.8 Mq, our first observational constraint from a half 
dozen white dwarfs in this single star cluster suggests 
the actual remnant masses are just slightly larger than 
the minimum core mass (with a slight extrapolation). 
This implies that the progenitors of such low mass stars 
had just enough mass to experience thermal pulses on 
the asymptotic giant branch and that the cores of these 
stars will not grow appreciably during this evolution. As 
M4 is a typical population II system in terms of its old 
age and low metallicity, it is reasonable to assume that 
these trends may extend to stars that are presently in 
similar phases of stellar evolution both in the general 
field halo population and in other globular clusters with 
similar properties to M4. 

5.2. The Dependence of Mass Loss on Metallicity 

In addition to extending the initial-final mass map- 
ping to lower masses, M4 represents the first star cluster 
with a significantly sub-solar metallicity in which white 
dwarf masses have been measured. The mean metallicity 
of the cluster is [Fe/H] = -1.10 (Marino et al. 2008), 13 
times more metal-poor than the bulk of the open clusters 
used to define the relation and 40 times more metal-poor 
than the metal-rich system NGC 6791 ([Fe/H] = -^0.40, 
Peterson & Green 1998). For this cluster, after account- 
ing for the new Stark broadening calculations discussed 
earlier, Kalirai et al. (2007) actually measured a mean 
mass for nine white dwarfs of Mgnai = 0.46 ± 0.06 Mq, 
much lower than expected. As they show, the bulk of 
these white dwarfs contain helium-cores, after evolving 
through a channel involving intense mass loss on the red 
giant branch, likely driven by the high metallicity of the 
cluster. The single data point shown in Figure 9 is based 
on the only clearly non He-core white dwarf in that sam- 
ple, and does not affect the discussion above. However, 
the results from NGC 6791 suggest strong dependen- 
cies of mass loss on metallicity, especially at super-solar 
metallicities. 

To date, we lack a good theoretical understanding 
of the dependence of post main-sequence mass loss on 
metallicity. The general expectation in evolutionary 

models is for more efficient transport of material outward 
from the star if the metallicity is higher, and therefore 
higher rates of mass loss. However, the exact treatment 
depends on the modeling of winds on the red giant branch 
and thermal pulses on the asymptotic giant branch (e.g., 
see Weidemann 2000 and Habing 1996), phases of stellar 
evolution that are poorly understood. As Catelan (2008) 
summarizes, mass loss prescriptions along the red giant 
branch alone show very different behaviors with metal- 
licity. For example, over a 1 dex metallicity change from 
[Fe/H] = -1.0 to 0.0, the MuUan (1978) and Schroder & 
Cuntz (2005) laws predict that AA/rcjb should increase 
by <0.05 Mq whereas the Goldberg (1979) and Judge 
& Stencel (1991) laws suggest AM increases by 0.15 - 
0.25 Mq. 

To qualitatively compare the dependence of mass loss 
in post main-sequence evolution with metallicity to these 
estimates, we can first parameterize the initial-final mass 
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relation. The solid line in Figure 9 is a linear fit to the en- 
tire sample of open cluster white dwarf - main-sequence 
progenitor mass pairs from Kalirai et al. (2008), 

Mfinai = (0.109 ± 0.007) Minitiai + 0.428 ± 0.025 Mg. 

The zero point of this relation has been adjusted to 
reflect a global offset of +0.034 Mq in the white dwarf 

mass measurements, although the slope may in fact also 
be affected if the offset is found to vary with mass (see 
Figure 11 and 12 in Tremblay & Bergeron 2009). The 
bounds on this relation range from the oldest open cluster 
white dwarf in NGC 6791 to the youngest white dwarfs in 
the Pleiades, and therefore are Minitiai = 1-1 ^ 6.5 Mq. A 
small extrapolation of the empirical relation to the mass 
of the present day turnoff in M4 (e.g., ^0.8 Mq) yields 
Affinal = 0.52 ± 0.02 Mq, which is therefore consistent 
(at l(j) with our measured value of the masses of M4's 
white dwarf population. These data therefore hint that 
mass loss rates have a weak or no dependence on mctal- 
licity, at least over the range extending from metal-poor 
populations such as M4 to roughly solar metallicity. In- 
terestingly, we know from our study of NGC 6791 that 
mass loss rates are strongly correlated with metallicity 
at extremely high [Fe/H], and therefore this may lead to 
an interesting shape for the empirical mass loss metal- 
licity relation from white dwarf spectroscopy. Of course, 
stronger conclusions will benefit greatly from additional 
observations of clusters spanning a wide range of metal- 
licities, especially at the metal-poor end. For comple- 
tion, if we add our new data point for M4 to the initial- 
final mass relation and recalculate the best fit ignoring 
any possible metallicity-related biasses, wc find a slightly 
shallower relation with Mgnai = (0.101 ± 0.006) MinUiai 
+ 0.463 ± 0.018 Mq. 

Finally, we note that the metallicity of M4, [Fc/H] = 
— 1.10, is only slightly more metal-rich than the global 
metallicity of both the stellar halo of the Milky Way 
([Fe/H] = -1.6, a = 0.6; Morrison et al. 2003) and M31 
([Fe/H] = -1.5, cr = 0.7; Kalirai et al. 2006). Assuming 
the dominant population in these galactic components is 
old (e.g., >10 Gyr), the same 34% steUar mass loss frac- 
tion should be representative of stars evolving in these 
populations. By directly inputing this mass loss into gen- 
eral stellar evolution models of similar age and metallicity 
(e.g., Girardi et al. 2000), we can add an important ob- 
servational constraint to the modeling of light from the 
stellar halos of spiral galaxies using population synthesis 
techniques (e.g., Bruzual & Chariot 2003). 

6. THE AGES AND DISTANCES OF GLOBULAR CLUSTERS 

In Section 1 we introduced the extensive HST observa- 
tions that have defined some of the deepest CMDs ever 
measured, with any instrument (e.g., the HST Cycle 13 
observations of NGC 6397 in Richer et al. 2008). Al- 
though these data sets have led to a wide range of astro- 
physical studies, the primary goal of the observations has 
aimed to accurately establish turnoff-independcnt white 
dwarf cooling ages for population II star clusters. In this 
respect, the present study serves to establish several of 
the input ingredients that are needed to make this mea- 
surement. For example, knowledge of the masses and 
spectral types of white dwarfs, and the initial-final mass 



relation, are all required to fully model an observed white 
dwarf cooling sequence (e.g., Hansen et al. 2004; 2007). 

6.1. Specific Ingredients for Accurate Age Measurements 

For M4, two recent studies have modeled the white 
dwarf cooling sequence using independent data sets, and 
both have concluded with similar ages (however, see dis- 
cussion below). First, Hansen et al. (2004) used a full 
2-D treatment of the white dwarf cooling sequence in 
the observed HST WFPC2 CMD (i.e., luminosity and 
color information are used), and just recently, Bedin et al. 
(2009) have used a more simplified technique of age dat- 
ing based solely on the luminosity function in newer and 
deeper ACS data. An analysis of these studies illumi- 
nates several interesting choices relative to the findings 
in this paper: 

• The present day mass of cluster white dwarfs form- 
ing today sets the normalization for the masses 
along the entire cooling sequence. In their mod- 
eling, Hansen et al. (2004) assume the mass at the 
tip to be 0.55 Mq and Bedin et al. (2009) assume 
0.54 Mq. As we have shown, the actual mass of 
white dwarfs at the bright end of the cooling se- 
quence is only slightly lower at 0.53 Mq, and there- 
fore the mass distribution (and hence cooling rates) 
along the entire cooling sequence are not very dif- 
ferent from the assumptions in these studies. 

• Hydrogen (DA) and helium (DB) atmosphere white 
dwarfs have different cooling rates, and so the frac- 
tion of stars cooling with these spectral types leads 
to structure in the luminosity and color function 
that must be modeled correctly to infer the age 
of the population from white dwarf cooling theory. 
For M4, NGC 6397, and NGC 6752, spectroscopy 
of the white dwarf cooling sequence indicates that 
100% of the stars at the bright end have hydrogen 
atmospheres (see Moehler et al. 2004 and Davis 
et al. 2009). Hansen et al. (2004) treated the frac- 
tion of DB white dwarfs as a free parameter and 
foimd the best fit to yield an upper limit of 40% 
(good fits were obtained as long as more than 60% 
of the white dwarfs are DA), whereas Bedin et al. 
(2009) choose an arbitrary value of 30% (e.g., the 
field disk ratio). 

• The initial-final mass relation is required to prop- 
erly evolve the progenitor stars on to the white 
dwarf cooling sequence. The normalization and 
form of the initial-final mass relation adopted in 
both of these studies (e.g.. Wood 1992), is based 
on empirical constraints from solar metallicity data 
points at higher masses, and a theoretical anchor 
at low masses. This can now be updated to our 
new result given earlier. 

We encourage future theoretical modeling efforts to 
make hill use of this new information in their analysis, 
preferably on the newer ACS imaging data of M4. Of 
course, the white dwarf models themselves should also 
be updated with new physical insights. For example, the 
contribution to the opacity of cool white dwarfs from the 
far red wing of the Lya line (Kowalski 2007) alone can 
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yield ages of globular clusters (e.g., NGC 6397) that are 
higher by -0.5 Gyr (Kowalski 2007; B. Hansen 2009, 
private communication) . 

6.2. A Closer Look at M4 's White Dwarf Cooling Age: 
Model and Distance Uncertainties 

We noted above that the Hansen et al. (2004) and Be- 
din et al. (2009) studies concluded with similar white 
dwarf cooling ages for M4. Formally, Hansen et al. (2004) 
reported an age of 12.1 Gyr {2a lower limit of 10.3 Gyr) 
and Bedin et al. (2009) measured an age of 11.6 ± 0.6 Gyr 
(internal errors). A deeper look at these two studies sug- 
gests they are, in fact, grossly inconsistent with one an- 
other given differences in white dwarf cooling models and 
uncertainties in the cluster distance modulus. As our 
present work is closely linked to obtaining very accurate 
absolute age measurements of globular clusters (through 
the study of white dwarfs), we briefly summarize these 
two studies below and highlight the respective assump- 
tions and differences that are being made. 

In the Hansen et al. (2004) study, the authors compare 
their data to three sets of models and adopt a distance 
to M4 based on the well established technique of sub- 
dwarf distance fitting. They find d = 1.73 ± 0.10 kpc 
for M4 (e.g.. Richer et al. 1997), which is identical to 
independent distance measurements for the cluster from 
the Baade-Wesselink distance using the M4 RR Lyrae 
stars {d = 1.73 ± 0.01 kpc, Liu & Janes 1990) and 
from astrometry (d = 1.72 ± 0.14 kpc, Peterson et al. 
1995). At present, it appears the latter estimate is the 
most secure. For the subdwarf measurement, the com- 
parison was based on pre-Hipparcos parallaxes, and for 
some other clusters, a calculation with post-Hipparcos 
parallaxes indicates a longer distance by ~0.2 mags than 
initially estimated (Gratton et al. 1997). Similarly, the 
Baade-Wesselink method yields luminosities for M4's RR 
Lyrae stars that are fainter by ~0.2 mags as compared to 
the My-[Fe/H] relation derived by Gratton et al. (1997) 
for clusters with low extinction (using post-Hipparcos 
main-sequence fitting techniques). Unfortunately, re- 
deriving the distance to M4 using any photometric tech- 
nique is problematic given the very high extinction along 
the line of sight {Ay > 1), the variation in Ay across the 
face of the cluster, and the non-standard shape of the 
extinction curve along this line of sight (which passes 
through the Scorpius-Ophiucus dark cloud complex - 
Vrba, Coyne, & Tapia 1993). Irrespective of the adopted 
distance (although see below), Hansen et al. (2004) find 
that the white dwarf cooling age they derive for M4 from 
their own models (e.g., Hansen 1999) is consistent with 
that measured using the Chabrier et al. (2000) models, 
however they rule out the Salaris et al. (2000) models 
which yield ages that are older than the Universe. There- 
fore, this suggests that, for the same distance modulus 
and data set, simply changing the white dwarf model 
yields an absolute age difference of 2 Gyr. 

In the recent Bedin et al. (2009) study, the authors 
use only the Salaris et al. (2000) white dwarf models, 
and still conclude with an age similar to that derived by 
Hansen et al. (2004). The difference is actually reflected 
in the distance to M4; Bedin et al. (2009) use a much 
larger (i.e., 15%) value of d = 1.98 kpc. This value is 
formally inconsistent with the three independent studies 
above. Bedin et al. (2009) measure the distance to M4 



by simply fitting a single isochrone of approximately the 
correct metallicity to the turnoff and horizontal branch. 
In principle, this method can of course yield the distance 
to any stellar population with a measured turnoff, how- 
ever it has almost always been avoided for a large number 
of reasons. First, and most importantly, there are uncer- 
tainties in both the microphysics built into the stellar 
evolution models (e.g., opacities, equation of state ef- 
fects, and nuclear reaction rates) and in the treatment of 
processes that do not come from first principles (e.g., con- 
vection, core-rotation, diffusion, gravitational settling, 
etc.).^'* Because of these unknowns, simply adopting an 
isochrone from another groups' models would lead to a 
different measured distance modulus for the cluster (e.g., 
different models predict different horizontal branch lumi- 
nosities). Second, such comparisons of data and models 
to yield distances require transformation of the model 
from the theoretical to the observational plane, result- 
ing in another source of error. Finally, in this specific 
case, this comparison does not alleviate any of the afore- 
mentioned uncertainties related to our lack of knowledge 
of the reddening and extinction along this line of sight. 
In the end, the application of the Salaris et al. (2000) 
models and the longer distance modulus derived by Be- 
din et al. (2009) effectively offset each other, yielding a 
derived age similar to that of Hansen et al. (2004). 

Summarizing, there are clearly both intrinsic differ- 
ences in white dwarf models calculated by different 
groups and there is a lack of consensus on the correct 
distance to M4. For the former, increased testing of the 
white dwarfs models on star clusters with well known 
fundamental parameters (e.g., distance, metallicity, and 
foreground reddening) are sure to resolve the correct pre- 
scriptions. For the latter, it appears that only the Peter- 
son, Rees, & Cudworth (1995) astrometric measurement 
is free of the unknown reddening and extinction laws 
along this line of sight, and therefore should be preferred. 
Until these differences are unambiguously resolved (e.g., 
the Space Interferometry Mission could directly measure 
the distance to M4 to within a few percent), the error 
budget from these uncertainties will remain large and the 
new physical constraints on masses, spectral types, and 
the initial-final mass relation are secondary to the analy- 
sis. Of course, if the properties of bright white dwarfs in 
other population II systems are similar to M4, the new 
results can have immediate significance on those studies. 
This is especially exciting given the remarkable data set 
that has been presented for NGC 6397 by Richer et al. 
(2008) and will soon be available for 47 Tuc (H. Richer, 
Cycle 17 ACS/WFC3 Proposal GO-11677). 

6.3. White Dwarfs as Standard Candles 

The mass measurement of white dwarfs in a globular 
cluster provides important leverage to the use of these 

stars as distance calibrators. Although the white dwarf 
cooling sequence is fainter than the bright part of the 
main-sequence in a star cluster, uncovering the sequence 
can offer several advantages over canonical techniques 
involving main-sequence distance fitting. For example, 
the location of the bulk of the white dwarf cooling se- 
quence in a CMD does not depend on the metallicity or 

VandenBerg, Stetson, & Bolte (1996) present a detailed dis- 
cussion of these effects on stellar evolution models. 
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age of the population, or on any of the other systemat- 
ies that plague main-sequence evolution (e.g., convection 
theory). However, the location does depend on mass 
(e.g., Wood 1995), with a 0.1 Mq difference translating 
to a A(mag) = 0.25 shift (i.e., a 12% distance offset). So, 
first, a simple offset can now be calculated to correct the 
location of any local calibrating sample of white dwarfs 
with parallax measurements (usually 0.60 Mq for field 
white dwarfs) to the expected location of 0.53 Mq stars. 
Following this, a straightforward comparison with the 
apparent luminosity of the globular cluster white dwarf 
cooling sequence can yield the distance to the popula- 
tion (e.g., see Renzini et al. 1996). Unfortunately, it is 
difficult to actually make this measurement and inde- 
pendently derive the distance of M4 given the aforemen- 
tioned uncertainties in the slope of the extinction curve 
along this line of sight. However, the true distances of 
other clusters with well constrained reddenings and ex- 
tinctions can be easily calculated from their cooling se- 
quences and the new mass measurements. 

7. SUMMARY 

We have presented the first directly determined spec- 
troscopic mass measurements of white dwarfs that are 
known to be population II members. Our measurements 
suggest that the final remnant mass at the tip of the 
white dwarf cooling sequence in the globular cluster M4 
is Mfinai = 0.53 ± 0.01 Mq, in nice agreement with expec- 
tations based on stellar evolution theory. These results 
suggest that low mass population II stars will lose ~34% 
of their mass through post main-sequence stellar evolu- 
tion. If, in fact, all Milky Way globular clusters have 
the same age, then this mass loss rate should be appli- 
cable to the evolution of stars in all but the most metal- 
rich Galactic globular clusters. We use this measurement 
to first extend the initial-final mass relation of stars to 
the lowest possible limit, given the age of the Universe. 
This single observation also hints that mass loss rates are 
not significantly different for metal poor stars when com- 
pared to solar metallicity stars of the same mass. Finally, 
we offer several pieces of new input that are required to 
measure more accurate white dwarf cooling ages of glob- 
ular clusters in the Galactic halo. 

Imaging observations of globular clusters with ground 
and space based telescopes have now uncovered thou- 
sands of white dwarfs in population II clusters (e.g., see 
Figure 1). The science that results from these purely 
photometric surveys is greatly enhanced with spectro- 
scopic follow up, which we have shown can yield the 
spectral types of the remnants (Davis et al. 2009), the 
mapping of initial to final mass, and the mass at the tip 
of the observed cooling sequence (and therefore a cali- 
bration for the masses along the entire cooling sequence) . 
All of the nearest globular clusters are in the southern 
hemisphere, and given their larger distances compared to 
M4, an extension of this work to these systems will re- 
quire superior blue sensitivity and throughput, and much 
larger allocations of spectroscopic time on 8-meter class 
telescopes (e.g., see Moehlcr et al. 2004). Future white 
dwarf spectroscopy in these clusters will benefit greatly 
from a 30-m telescope with a multiobject spectrograph 
with the bluest sensitivity. Such a telescope would not 
only allow a direct measurement of the mass loss rates 
in an expanded set of open and globular clusters, but 



would permit the first initial-final mass relations to be 
built from observations of a set of white dwarfs in a 
single cluster. Deep spectroscopy of white dwarfs over 
a range in luminosity (and therefore mass) can be en- 
gaged in a given cluster, minimizing the systematics in- 
troduced by metallicity variations in samples of clusters. 
Fortunately, the nearest three globular clusters also sam- 
ple a wide range in metallicity (e.g., a factor of 20 from 
NGC 6397 at [Fe/H] ^ -2.0 to 47 Tuc at [Fe/H] = -0.7), 
and therefore the dependence of mass loss on metallicity 
can be probed in detail with such observations. 
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Fig. A1. — These snapshot figures show several pieces of information for each of the white dwarf candidates WD 00 - WD 08. The two 
panels on the left illustrate the location of these stars on the ground based CMD (top) and the HST CMD (bottom), where available. The 
individual candidate star is marked in each panel with a larger red diamond. The location of the star relative to our imaging field of view 
is also shown in the top-right panel (see Figure 2 for details on this plot). In the bottom-right panel, we present a small snapshot image of 
the candidate. Several objects that are targeted clearly have very little flux in the ground based image and/or are located close to bright 
stars. These stars are HST detections and are discussed, and removed from subsequent analysis, in § 3. Note, the identifications of "LRIS" 
in this figure correspond to "WD" in Table 1. A higher resolution version of this figure is available in ApJ. 
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Fig. A2. — Same as Figure Al for WD 09 — WD 17. A higher resolution version of this figure is available in ApJ. 
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Fig. A3. — Same as Figure Al for WD 18 — WD 26. A higher resolution version of this figure is available in ApJ. 
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Fig. A4. — Same as Figure Al for WD 27 - WD 30, and Gem WD 05, 11, 15, and 17. A higher resolution version of this figure 
available in ApJ. 
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